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Abstract 

Saltation is an important geological process and the primary source of atmospheric mineral dust aerosols. Unfortunately, no 
studies to date have been able to precisely reproduce the saltation process because of the simplified theoretical models 
used. For example, sand particles in most of the existing wind sand movement models are considered to be spherical, the 
effects of the sand shape on the structure of the wind sand flow are rarely studied, and the effect of mid-air collision is 
usually neglected. In fact, sand grains are rarely round in natural environments. In this paper, we first analyzed the drag 
coefficients, drag forces, and starting friction wind speeds of sand grains with different shapes in the saltation process, then 
established a sand saltation model that considers the coupling effect between wind and the sand grains, the effect of the 
mid-air collision of sand grains, and the effect of the sand grain shape. Based on this model, the saltation process and sand 
transport rate of non-spherical sand particles were simulated. The results show that the sand shape has a significant impact 
on the saltation process; for the same wind speed, the sand transport rates varied for different shapes of sand grains by as 
much as several-fold. Therefore, sand shape is one of the important factors affecting wind-sand movement. 
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Introduction 

Saltation plays a key role in various geological processes, 
including wind erosion, sediment transport, and the formation of 
sand dunes [1]. In-depth studies of sand saltation, which is the 
main manifestation of sand movement near the bed surface 
[2], [3], are considered to be essential to comprehensively 
understanding wind erosion and desertification processes and, 
furthermore, to preventing sand disasters [4] . Many scholars have 
carried out studies to simulate the saltation process [5-14]. 
However, significant discrepancies still exist between the simulated 
results and measurements [1]. Current models of saltation make 
many simplifications because of the complexity of sand saltation. 
For example, most models assume that sand particles are 
spherical, and the effect of mid-air collision is ignored. 

In fact, spherical and sub-spherical sand grains generally 
comprise no more than 30% of the total sand [15]. Many scholars 
have studied sand shape in the desert [16-20]. For example, Sagga 
[16] and Juan [17] studied the roundness of sand grains from 
dunes in Saudi Arabia and Mexico, respectively. Sagga found that 
the sands of the inter-dunes are more rounded than are those of 
the adjacent dunes and that sand particles tend to be less rounded 
with increasing dune height. Juan proposed that aeolian and 
marine processes might generate different grain sizes with dilferent 
compositions. Cheng et al. analyzed aeolian sandy physical 
characteristics and measured the roundness, average particle size 
and size distribution of desert sand north of Beijing. They 



indicated that these parameters have an important influence on 
the interactional force of aeolian sand [18]. 

Researchers have investigated the mechanical properties of non- 
spherical particles for over half century. Barton studied the slip 
correction factors for non-spherical bodies in continuum flow in 
1973 [21]. Then, Gavze and Shapiro investigated particles in a 
shear flow near a solid wall considering the impact of non- 
spherical grains [22]. In recent years, more studies have been 
conducted in this field. For example, Yin et al. developed a model 
of motion of cylindrical particles in a non-uniform flow [23]. Xu 
et al. discussed drag and lift forces of rotational non-spherical 
particles [24]. Zastawny et al. considered the drag and lift force 
and torque coefficients for non-spherical grains in flows [25]. 

Based on the above studies of the motion of particles in a shear 
flow, our study first proposes a relationship between sand shape 
and the threshold friction velocity of wind. Then, a sand saltation 
model, which takes into consideration the coupling effect between 
wind and sand grains, the effect of the mid-air collision of sand 
grains, and the effect of the shape of sand grains, is established, 
and the saltation process of non-spherical particles is simulated. 
Finally, the effects of different shapes of non-spherical sand grains 
on wind-sand movement are analyzed in detail in this paper. 

Methods 

1. Analysis of sand shapes 

Blott and Pye [26] applied Beckman-Coulter to analyze natural 
sand particles, including desert dune and coastal dune sand 
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particles, and found that the average roundness *F of natural sand 
is in the range of 0.82-0.86. The average roundness I* is defined 
as the ratio of the surface area of a sphere to the surface area 
Ap of a sand particle with an equivalent volume, that is 

^ = Ay/A„ = nd';/Ap (1) 

where rf,, is the equivalent diameter of the sphere with the same 
volume as that of the sand particle. 

Let 5,(a,b,c) represent non-spherical particles of different 
shapes, where a, b and c are their size parameters, respectively. 
In this study, we examined the effects of four types of non- 
spherical particles, including ellipsoid, cube, cylinder and frustum, 
on wind-sand movement. 

1) Ellipsoid-shaped sand particles: particles with a shape 
approximate to an ellipsoid are simplified to an ellipsoid, 
with a, b and c indicating the lengths of the long axis, mid 
axis, and short axis, respectively. We selected two types of 
spheroids, S\{l,l,ll / 5) and 82(1, 1/3,1/ 3), respectively. 

2) Cube-shaped sand particles: particles shaped like cubes are 
reduced to a regular cube SjiUlJ). 

3) Cylinder-shaped sand particles: particles shaped like cylinders 
are simplified to a regular cylinder, with a and b indicating the 
diameter of the bottom and c indicating the height. Two 
cuboids were selected, Sniil/ 5,31/ 5,1) and 8^(1,1,31 / 5); that 
is, their bottom diameters are 31/5 and /, respectively, while 
their heights are / and 31/5, respectively. 

4) Frustum-shaped sand particles: irregular frustum-shaped sand 
particles are reduced to a regular frustum, wherein a, b and c 
indicate the upper end diameter, the lower end diameter, and 
the height, respectively. The regular frustum Sf,{l/l,l,l) was 
selected with a, b and c corresponding to 1/2, I, and /, 
respectively. 

The symbol / in the above 1) to 4) is a constant. To consider the 
effects of different shapes, three shape parameters are introduced 
[27] as follows: 

1. The aspect ratio, £ = <i|| /rfi , where d\\ and c/i are the 
maximum diameters or sizes of the horizontal and vertical 
directions of the non-spherical particles (Table 1); 

2. The ratio of the surface areas, A*^^^.j- = Asurf/(^df); Asmf is the 
surface area of non-spherical grains. By comparison with 
formula (1), ^*,^^. = 1/>I' (Table 1); 

3. The ratio of the projected areas, A*i„-oj = '^^proi / ^d^', ^pwj is the 
projected area of the non-spherical grains (Table 1). 



2. Drag coefficient of non-spherical grains 

The flow field can be divided into three parts according to the 
Reynolds number [28]: a creeping field of flow (Re«l), an 
intermediate field of flow and a Newton field of flow (ca. 
2000 < Re < 300000). 

Under the actions of the three fields of flow, the drag coefficient 
of spherical particles often is [2], 

Cz, = 24/Re + 6/(1 +(Re)'''2) + 0.4 (2) 

where the Reynolds number 'R.e = DpVr/ v. 

To consider the impact of particle shape, the modified 
coefficient of shape is defined as [27], 

where Co.shape is the drag coefficient of a non-spherical sand 
particle, and CD,spliere is the drag coefficient of a spherical sand 
particle. 

Under different fields of flow, the modified coefficient of particle 
shape will be different. The following is an analysis of the shape- 
modified coefficients of non-spherical particles subject to three 
types of fields of flow: 

The particles in the creeping flow field [27], Re«\: 

CLpe = + 2/4^/3 (4) 

The particles in the Newton flow fleld [27], 
2000 <i?,< 300000: 

CL,e = 1 + 1 -^iAlrf - 1)'^' + 6.7(^L,, - 1)'^ ^ 1 (5) 



CLp. = 1 +0.7(^,V - 1)'^' + 2.4«,/ - 1),£> 1 (6) 

The particles in the intermediate flow field [27]: 
Q = Cz)/C,/,„,,,, Kp = ClhapeRepl Cl^ape' whcrc Q and Rip are 
the substitute of the intermediate flow field. Some scholars have 
obtained the following empirical formula [27]: 



Table 1. Parameters of sand particles with different shapes. 
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where S] and S, are two different ellipsoid-shaped sand particles, ^3 represents cube-shaped sand particles, and represent two types of cylinder-shaped sand 
particles, and S(, represents frustum-shaped sand particles. 
doi:1 0.1 371/journal.pone.01 05208.t001 
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Figure 1 . Schematic of analysis of the forces acting on a sand particle. Particles B and C are the supporting sand particles on the bed surface. 
doi:1 0.1 371 /journal.pone.01 05208.g001 



For sand particles with a circular cross-section: 
C* = 24[1 + 0. 1 5(R*^f^^'^] I Rl^, + 0.42/[l +42500 / (i^:^''"'] (7) 

For sand particles with a non-circular cross-section: 



Q =24[1 + QmiiKf] R% + QA1 [\ + 33 UrJ''] 



(8) 



3. Starting friction speed of non-spherical particles 

Before starting to move, sand particles on the bed surface are 
mainly subject to the gravity W, the drag force Fr,, the force of 
supporting the bed Fj^, and the cohesive force F,, among the sand 
grains, among which the gravity and the drag force are given as 
follows [11]: 



Fig. 1 shows a schematic of the force analysis among three non- 
spherical particles on the bed surface, in which particles B and C 
are the supporting sand particles of the bed surface. Relative to 
Point 0\, the sand particle A is subject to the moment of force. 
According to the moment balance principle, as the sand particles 
satisfy formula (11), they wiU leave the bed surface: 



FDLD>WLf. 



(11) 



where Lw = l cos d and Lb = / sin 9. 

From inequality (10), the critical condition for the sand particles 
to start saltation is: 



Fd smd=W cos d 



(12) 



Substituting formulas (9) and (10) into (12), we obtain the 
starting speed of flow as follows: 



W = Kd,'(pp-p)g/6 



(9) 



Uc = J Tid, (pp - p)g/ (3CdpA tan ( 



(13) 



FD = CDpAVj/l (10) 

where A is the maximum area of the particle cross-section. A wiU 
adopt different values for different sand shapes; Vr is the relative 
speed between the sand grains and the gas flow; p is the 
aerodynamic coefficient of viscosity; pp is the density of the sand 
particles; p is the density of air; and g is the acceleration due to 
gravity. 

Here, we focus on dry sand particles. Thus, for large particles, 
the cohesive force among the particles is very small and can be 
ignored [2]; that is, Z',, =0.0. 



For a given sand particle shape, is only related to 9, which is a 
random variable, depending on the acting point of the supporting 
particle against the saltation particle. 

4. Wind field equation 

Before sand particles start to move near the ground's surface, 
the motion of air flow can be described by the Navier-Stokes 
equation: 

pdu/dt + pu-W= -Vp + V-T + pg (14) 
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where u is the velocity of wind, t is time, g is the gravitational 
acceleration, f is the Reynolds' tangential stress, and p is the 
pressure. 

When sand particles leave the bed and move into the saltation 
layer, the blocking effect of the particles against the wind field will 
change their distribution, and this blocking effect can be 
equivalent to adding a correction term, the volume force, to the 
right of the Navier-Stokes equations: 



pdu/dt + pu-W= -Vp + S/-T + pg + F^ 



(15) 



Letting the direction of air flow movement be in the x-axis and 
the direction perpendicular to the bottom be the y-axis, m-Vm = 0, 
d/8y»8/dx, d/dy»d/dz. Then, we obtain the motion differen- 
tial equation of the wind field component in the x direction of in 
the saltation layer of the wind sand flow. 



p8u/8t=dT:/8y+Fj^(y) 



(16) 



If a two dimensional steady state is fully developed in the flow 
field, then the left term in Eq. (16) is equal to zero, that is, ^ =0; 
thus, Eq. (16) becomes 



0 = dT/dy+F^{y) 



(17) 



The blocking effect of saltating particles against the wind field 
changes with height, and thus the stress acting in the field of the 
wind also changes with height. Integrating Eq. (17) with respect to 
height results in Eq. (18): 



ymax 

^iy)+ I F^{y)dy = Xb 



(18) 



wind sand flow [13]: 
du/dy = 



Fxiy)dy / p 



, 1 

2 



ky (22) 



In accordance with the logarithmic distribution profile that the 
wind field obeys 



M = 5.75m, IgCf/jo) 



(23) 



where yQ = dy/7>Q is the roughness of the bed surface. 

The critical starting height is yc = d/2. Bringing u = Uc and 
y=yc into formula (23), we can obtain the critical starting friction 
speed Utt of non-spherical sand particles. 

5. Saltation equation 

The analysis of sand grains subject to forces indicates that in the 
process of their movement, they experience only the gravity W and 
the drag force Fo, as shown in Formulas (9) and (10), respectively. 
The basic equations of sand saltation movement are: 



pd^x/df=f],-- 



- pACd{x-u)\I {k-uf +f / 2 (24a) 



nipd^y/dt^ =fy - rripg =-pA CDy\J{x-uf+f j 1- ntpg (24b) 

in which rrtp is the mass of the sand particles, g is the gravitational 
acceleration, u is the wind speed in the x-direction, and x and y 
are the velocity components of the particles in the x-direction and 
y-direction, respectively. 

The starting initial velocity distribution of saltating sand 
particles was obtained by Dong et al. in their experiments [29]: 



where j'max is the max height to which the particle can move 
upward, is the shear stress of the air flow from the blocking 
effect from the sand particle above j'max, Tj, is the air flow shear 

ymax 

stress at height y, and J Fx(j)dy is the particle-loaded 

y 

tangential stress [7], [11]. 

Based on the Prandti mixing length, the Reynolds tangential 
stress becomes 



■z = p,du/dy 



(19) 



/("/.) =/("o,vo) =/(«o)/(i'o) 



(25) 



where /, A, w, Uc, fi, xq, Ai, t\, A2, and t2 are recession 
parameters. 

6. Resistance of saltating sand groups to the wind field 

The interaction between the wind flow and the saltation layer 
plays an important role in the saltation process [3]. After 
considering the sand shape, the resistance from sand particles 
experienced by the wind per volume can be described [13]: 



In the saltation layer of the wind-sand flow, the use of the 
equivalent friction speed can obtain [7], [11]: 

Uteff = \/tJp (20) 



p,=pkUteff (21) 

where u„ff is the equivalent friction speed in the saltation layer 
and varies with y, and Hi is the viscosity coefficient of the turbulent 
flow. Bringing Eqs. (19)-(21) into Eq. (18), we obtain the two- 
dimensional steady state flow field equation fully developed for 



F,(y) = s{nd,'' Pp/6)Y{vp){a,y(y) /\y^\+a,^{y) /\y^\)dvp (26) 

where a^^iy) and ax[(y) denote the accelerations of the particles 
moving upwards and downwards at height y, respectively, and y^^ 
and are the average speeds of the particles moving upwards and 
downwards at height y, respectively. 

7. Calculations of the sand transport rate (STR) and the 
SIR per width 

Considering the mechanism of sand collision, the STRs of 
irregular particles can be divided into two STRs representing 
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Figure 2. The drag coefficient of sand particles as a function of the Reynolds number. The relationship of the drag coefficient of sand 
grains of different shape to the Reynolds number is shown. Here, the black solid line indicates spherical sand grains, S\ and 5*2 are two different 
ellipsoid-shaped sand particles, 5*3 represents cube-shaped sand particles, S4 and represent two types of cylinder-shaped sand particles, and 5*6 
represents frustum-shaped sand particles. 
doi:1 0.1 371 /journal.pone.01 05208.g002 



moving particles with collision and those without collision in air. 
The STR at height y considering the aerial collision mechanism of 
sand particles and the co-initial speed distribution of starting 
particles is described as [13], [30]: 

QA}')= I smpf(vp)[(l -p^(v,,,y))Xfivp)/yj(\'i,)-(l -pl(\^,„y))x^(vp)/y^(V|,)]dV|, 

(27) 

+ I .impf(rp)\Pl(rp,y)ixj(Vp)/y^(Vi,))^ -pi(r,„y)(xi{v,,) / y^(v,,))Jdv,, 



Integrating Eq. (26) with respect to height, we obtain the STR 
per width of irregular particles [13], [30]: 

Qi«c= I I .sm,/(v,)[(l-p,(Vp,y))i,(v,,)/>',(iv)-(l-pi(Vp,>'))ii(i>)/>'i(v,,)|rfv,, 

„ (28) 

+ 11 ."«;/( >',)[PT(i'pJ')(it('';<) / y;(^p))^ -Pl(rp,y)(xi{Vp) j y ^(v p)) \dv 



where p\{yp,y) is the probability of collision of rising particles at 

height y with an initial lifting-oflF velocity of v^, (x|(Vp) j y\iyp))p is 

the ratio of the horizontal velocity to the vertical velocity after the 
aerial collision of the rising particles at height v with an initial 



lifting-oflF velocity of Vp considering the aerial collision mechanism, 
is the probability of collision of falling particles at 
heightj'with a co-initial lifting-oflf velocity of Vy,, and 

(X|(v,,)y/j|(v^))y, is the ratio of the horizontal velocity to the 

vertical velocity after the aerial collision of falling particles at 
height 31 with a co-initial lifting off velocity Vp. 

8. Criterion of wind sand flow up to dynamic equilibrium 

The critical impact starting velocity of wind is approximately 
80% of the starting wind velocity of a critical fluid [2]. From the 
starting wind speed of a critical fluid Uc obtained from Eq. (13), the 
critical impact starting shear stress considering the sand shape is 

t„ = 0.8t,/=0.8pmJ (29) 

where T„ and x,f are the critical impact starting shear stress of the 
wind sand flow and the starting shear stress of a critical fluid [13]. 

9. Constants in the model 

For the convenience of comparison, the equivalent diameters of 
various shapes of sand particles are assumed to be the same. We 
select two groups of sand particles in this article, dv = 0.25 mm and 
0.35 mm; the sand density is Pp = 2650 kg/m^, the air density is 
p= 1.22 kg/m'\ and the dynamic viscosity coefficient of air is 
v=1.5xlO^'; and we assume that the initial wind speed profile 
obeys a logarithmic distribution. 
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Figure 3. The critical starting friction wind velocity for sand grains of different shape, (a). Comparison of the critical starting friction wind 
velocity of non-spherical particles with that of spherical particles (equivalent diameter rf,. = 0.25;;im). (b). Comparison of the critical starting friction 
wind velocity of non-spherical particles with that of spherical particles (equivalent diameter rf,. = 0.35m;;;). The black solid line indicates the spherical 
sand grains, Si and Si represent two different ellipsoid-shaped sand particles, S3 represents cube-shaped sand particles, S4 and S5 are two types of 
cylinder-shaped sand particles, and S^ represents frustum-shaped sand particles. 
doi:1 0.1 371/journal.pone.01 05208.g003 
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Figure 4. The drag force of different shaped particles with height. All the sand particles of different shapes take off with the same horizontal 
initial velocity of V\ and vertical initial velocity of V2. The black solid line indicates spherical sand grains. Si and S2 are two different ellipsoid-shaped 
sand particles, S3 represents cube-shaped sand particles, S4 and S5 are two types of cylinder-shaped sand particles, and S^ represents frustum- 
shaped sand particles. 
doi:1 0.1 371/journal.pone.01 05208.g004 
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10. Computational steps 

1) Substitute the parameters of non-spherical sand grains 
A*^gj and E into Eqs. (4)-(8) to calculate the corresponding 
shape parameters and calculate the resistance coefficients of 
different shape of particles Cj), then bring the coefficients into 
Eq. (10) to obtain the drag forces of different shapes of sand 
particles F]). 

2) Find the values of parameters Zti and m;„ and calculate the 
initial friction wind speed K», then obtain the initial wind 
speed distribution and denote them u^(y) or u'(y). 

3) Select the testing value of the sand liftoff rate from the sand 
bed s, calculate the co-initial speed distribution function of the 
saltation particles from formula (25), and denote the 
calculated results /'(v^) or f'{Vp). 

4) Substitute u^(y) or u'(y) into Eq. (23) to calculate the saltation 
trajectories of various sand groups and mark the obtained 
trajectories as x\t) or x'(t), y^(t) or y'{t), x\t) or x'{t) and 
y\t) ory(f). Put s,f(Vp) (or/'(Vp)), x\t) (or x'(0), and j^O 
(or y'(t)) into formula (26) to find the reaction force of the 
particles at different heights of wind Fx'{y). 

5) Substitute F^'iy) into Eq. (22) to find the wind speed profile 
modified by saltation particles u'^^(y), then find the shear 
stress of the wind acting on the bed surface t' + '. 

6) If t'"*"' is greater than (smaller than) t,,-, choose the greater 
(smaller) s and repeat Steps 3-6 until |t' + ' — t„ | < 0.001. 

7) Substitute the values of s, f'(vp), x'(t) and y'{t) calculated 
according to Step 6 and satisfying |t'+' — Th] <0.001 into 



formulas (27) and (28) to fmd the STR and the STR per width 
considering sand shape changes with height. 



Results and Discussion 

The relationships of the drag coefficients of six non-spherical 
sand granules to the Reynolds number are shown in Fig. 2. From 
the figure, it is obvious that the drag coefficients of non-spherical 
sand particles are greater than those of spherical granules and that 
the greater the Reynolds number R^p, the more obvious the 
difference between the non-spherical and spherical granules. 

Fig. 3 sho^vs the relationship between the critical starting 
friction wind speed of non-spherical sand flow and the Reynolds 
number. It can been seen from this figure, the critical starting 
friction wind speed given by the Bagnold is a fixed value when the 
particle size is invariant, while the critical starting friction wind 
speed in this paper changes with the particle Reynolds number. 
And the critical starting friction wind speeds for different shapes of 
sand particles are different for the same Reynolds number. For 
example, when the Reynolds number R^p = 1 84, the critical 
starting wind speed of spherical sand particles is = 0.176 m/.? 
(for D = 0.35 mm) (Fig. 3-b), while the critical starting wind speed 
of non-spherical particles S2 is M,, = 0.273 w/i (for D = 0.35 mm) 
(Fig. 3-b), representing an increase of 35.53% compared to the 
spherical case. Thus the effect of the sand shape on the critical 
starting friction wind speed is evident. This is mainly due to the 
drag coefficient changes with the particle Reynolds number and it 
is related to the shapes of the particles. 

Fig. 4 shows that when the friction velocity is 0.49 m/s and the 
equivalent diameter is 0.35 mm, the distributions of the drag force 
are affected by the shapes of the sand particles taking off at an 
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Figure 5. The sand transport rates of sand particles witli different shapes as a function of height, (a) and (b) show the sand transport 
rates of sand particles with the same equivalent diameter of <:/,, = 0.25mm as a function of height but at different friction wind velocities, (c), (d) and (e) 
represent the sand transport rates of sand particles with the same equivalent diameter of J, = 0.35mm as a function of height but at different friction 
wind velocities. Here, the black solid line indicates the spherical sand grains, S\ and S2 are two different ellipsoid-shaped sand particles, S3 represents 
cube-shaped sand particles, ^4 and Ss are two types of cylinder-shaped sand particles, and Sf, represents frustum-shaped sand particles. 
dol:l 0.1 371/journal.pone.01 05208.g005 
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Figure 6. The sand transport rate per width for sand particles of different shape, (a). Comparison of calculated sand transport rates per 
width for different shapes sand particles with equivalent diameters of df = 0.25mm with experimental data. (b). Comparison of calculated sand 
transport rates per width for sand particles of different shape with equivalent diameters of d„ = 035mm with experimental data. The black solid line 
indicates spherical sand grains, Si and 5*2 are two different ellipsoid-shaped sand particles, represents cube-shaped sand particles, 54 and are 
two types of cylinder-shaped sand particles, and S(, represents frustum-shaped sand particles. 
doi:l 0.1 371/journal.pone.01 05208.g006 



initial speed (such as at a horizontal initial velocity of 0.6 m/s and 
a vertical initial velocity of 0.735 m/s) along the elevation. It can 
be observed that during the saltation process, the changes in the 
drag forces of different shapes of sand particles are consistent, first 
increasing and then decreasing with increasing saltation height 
during their ascent, while during their descent, the drag forces first 
decrease and then increase with decreasing saltation height. Thus, 
the drag force for different shapes of sand particles at the same 
height can vary greatly. For example, the maximum drag force of 
spherical sand particles is 2.65x10 N, while that of non- 
spherical sand particles Si is 7.5xl0~^ N; the latter is 2.72 times 
larger than the former. 

It can also be observed from Fig. 4 that the drag forces applied 
to the non-spherical particles whose aspect ratios E are greater 
than 1 .0 are smaller than those applied to the spherical particles, 
while the drag forces applied to the non-spherical particles whose 
aspect ratios E are smaller than or equal to 1.0 are larger than 
those applied to the spherical particles. In addition, the smaller the 
aspect ratio, the greater the drag force applied to the sand grains 
and the smaller the maximum height the grains can reach after 
taking off. Because the drag force suffered by sand grains during 
their ascent has a vertical downward component, the maximum 
height they can reach is less than 0.0276 m (the maximum height 
the freely lifting motion at a vertical initial velocity of 0.735 m/s 
can achieve is 0.0276 m). The greater the drag force, the greater 
the vertical downward component and the smaller the maximum 
height the sand particles can achieve. For example, the drag forces 
acting upon non-spherical sand particles Si and S5 are smaller 
than those acting on spherical sand particles; the maximum height 
the spherical sand particles can reach is 2.15 cm, while the non- 
spherical particles Si can only reach a height of 1.45 cm. The 
shape of the sand grains therefore has significant impact on the 
drag force and the maximum height these particles can achieve. 
Thus, it also has a great impact on the structure of wind-blown 
sand. 

Fig. 5 shows the STR per unit area of different shapes of non- 
spherical particles with equivalent diameters of 0.25 mm and 
0.35 mm as a function of the saltation height. The solid line 
represents the distribution of the STR per unit area of the 
spherical sand particles as a function of height. As shown in Fig. 5, 



the variations in the STR per unit area for different shapes of sand 
particles with height are equivalent; that is, the STR per unit area 
decreases with increasing height. However, the STRs per unit area 
for different shapes of sand particles at the same height vary 
greatly. In general, the STRs per unit area of the sand particles ^2 
and 55 at various elevations are slightly larger than that of 
spherical sand particles of equivalent diameter. It can also be 
observed in Fig. 5 that the maximum heights of the wind sand flow 
formed by different shapes of sand particles at the same wind 
speed vary greatiy, even several-fold. For example, when the 
equivalent diameter is 0.25 mm and the friction velocity is 
0.49 m/s (Fig. 5-a), the sand transport rates per unit area for the 
non-spherical particles Si, S2, S3, S4, S^ and S(, at heights above 
13.2 cm, 34.8 cm, 27.6 cm, 30 cm, 37.2 cm and 20.4 cm, 
respectively, are equal to zero, while that of the spherical sand 
particles at a height above 34.8 cm is zero. Comparing the non- 
spherical sand particles Si and S5, the maximum height of S5 is 
2.82 times larger than the maximum height of Si. 

Fig. 6 shows the relationship of the STR per width of non- 
spherical sand particles to the friction wind velocity, among which 
Fig. 6-a and Fig. 6-b correspond to equivalent particle diameters 
of 0.25 mm and 0.35 mm, respectively. It is clear from the figures 
that the STRs per width for the different shapes of sand particles 
with friction wind speed exhibit similar variations; that is, they 
increase with increasing friction speed, in agreement with Dong's 
experimental results [31]. Carneiro et al. pointed out that the mid- 
air collisions exert an enormous effect on saltation flux and it 
depends strongly on restitution coefficient and wind speed [32]. As 
shown in Fig. 6, the impact of the sand shape on the STR per 
width is also significant. Under the same conditions, the STR per 
width for different shapes of sand particles can differ by several- 
fold. For example, when the equivalent particle size is 0.35 mm 
and the friction wind speed is 1.016 m/s (Fig. 6-b), the STRs per 
width of the non-spherical particles S5 and S(, are 0.056 kg /m/s 
and 0.015 kg/m/s, respectively, where the former is 3.73-fold 
larger than the latter. 

Thus, in the discussion of the STR, the sand shape is again an 
important factor. The main reason is that the air drag forces acting 
on different shapes of sand particles taking off at the same vefocity 
as well as their maximum heights are different, resulting in large 
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variations in energy when they fall to the bed surface. Thus, the 
differences in the STR are also significant. 

Conclusions 

In this article, we established a relationship between the starting 
friction wind velocity and the sand shape, simulated the STR and 
STR per width of wind sand flow reaching the steady state 
considering the sand shape and the distribution of co-velocities at 
which the sand particles lift off, and calculated the drag forces and 
the starting friction wind velocities of different shapes of non- 
spherical sand particles in wind sand flow. The results show that: 

(1) The drag coefficients of non-spherical sand particles are 
greater than those of spherical granules. The greater the Reynolds 
number is, the more obvious the difference between the non- 
spherical and spherical particles. 

(2) The change in the drag force is consistent for different shapes 
and sizes of sand particles; i.e., during their ascent, the drag force 
first increases and then decreases with increasing saltation height, 
while during their descent, it first decreases and then increases 
with decreasing saltation height. The drag force to which the sand 
particles are subjected is inversely proportional to the maximum 
height they can reach after their takeoff and closely related to the 
aspect ratio. 
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